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Abstract 

We discuss the phenomenological implications of assuming a Veneziano-type 
spectrum for the vector and axial-vector two-point functions in QCD at large 
N c . We also compare the phenomenological results with those of Lowest-Meson 
Dominance, and find that they agree. 



1 Introduction 



Quantum Chromodynamics simplifies when we take the number of colors, N c , to 
infinity ||T|; we will denote this limit by QCD^. Assuming confinement, it is generally 
believed that the large-iV c solution qualitatively resembles the real world enough to 
be taken as a serious starting point for a systematic approximation by expanding in 
inverse powers of N c 0. 

Unfortunately, no solution to QCD^ has been found to date, except in the two- 
dimensional case |J. We know that even in the large- N c limit this solution has to be 
very complex. Large N c , together with confinement and asymptotic freedom, predicts 
the existence of an infinite set of zero- width meson resonances in any channel with given 
quantum numbers. These resonances contribute to any Green's function made of quark 
bilinear operators corresponding to those quantum numbers. At very large energies a 
simplification occurs, as the infinite tower of resonances has to become equivalent to 
the free-quark loop of the parton model obtained in leading-order perturbative QCD 
(pQCD), with the connection provided by the operator product expansion (OPE). 
However, how this happens is not known, except in the two-dimensional case. We 
would therefore like to propose a simple model for this equivalence, which despite its 
simplicity is phenomenologically quite successful. 

In this work we will deal only with two-point functions in the chiral limit. We 
limit ourselves to this choice because these are the simplest Green's functions which 
obey non-trivial chiral symmetry and OPE constraints. The model is based on the 
observation that, although the equivalence between the resonance description of QCD^ 
and the quark-gluon description of pQCD is valid at asymptotically large energy scales, 
phenomenology suggests that asymptotic behavior may set in already at values of 
s = Q 2 ~ (1 — 2) GeV 2 . In the vector channel, the p meson is clearly lighter than 
this scale s, but all higher resonances (p',p",p" f , etc.) appear to lie already in this 
asymptotic region. Therefore, it is reasonable to consider this channel as consisting 
of a separate p resonance of mass M p , plus an infinite tower labeled by an index n, 
of masses My(n), which corresponds to the perturbative description in a way to be 
described below. In the axial-vector channel, on the contrary, the a\ may already be 
heavy enough to group all resonances together in a corresponding tower, M^(n). Here 
the exception is of course the Goldstone meson which, because of chiral-symmetry 
breakdown, is massless and stays away from the tower. 

How these resonances are spaced as a function of the index n is not known in 
QCDoo . Instead, we will rely on Regge theory [|J. It has been long suspected that 
there is a connection between large- N c and Regge theory 0, [| , but the precise link has 
never been found |J . One of the natural features of Regge theory is the existence of so- 
called daughter trajectories. Their most important feature for us is that they give rise 
to towers of resonances with a given spin obeying an equal-spacing rule with respect 
to the index n, i.e. the mass squared is given by M 2 (n) = a + bn. In two dimensions, 
where the solution to QCD^ is known, it can be shown that indeed M 2 (n) ~ n for 
large n. Here, in four dimensions, we will conjecture that this linear behavior in n is 
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also valid for both the vector and axial-vector towers, and assume that 



M^ A {n) = m 2 VA + n k 2 VA , (1.1) 

is a reasonable approximation for all mesons except the p meson and the pion. 

For simplicity, we will content ourselves with this equal-spacing (ES) ansatz pro- 
ducing only the leading term in pQCD, the parton-model logarithm. A possible remedy 
to this is to delay the onset of the equally spaced tower by including more and more in- 
dividual resonances which are not part of the tower, and thus move the "onset" of this 
tower to a higher energy, where the parton-model logarithm is a better approximation. 
In this way the description would become more accurate. Alternatively, one could take 
the n dependence in eq. ( |1 . 1|) to be more complicated!] such that non-leading pertur- 
bative corrections are reproduced from our ansatz as well. The price to pay, of course, 
is that this introduces more unknown parameters into our ansatz, which would have 
to be fixed by matching to the OPE and chiral perturbation theory (ChPT) to higher 
orders in a way similar to what we will describe for our simple ansatz in section 2. 

Much of this has been said before, although perhaps not within the unifying frame- 
work of the large- N c expansion. For instance, in the pioneering work by Bramon, Etim 
and Greco and Sakurai |7j , it was already observed that a resonance tower including 
as the first state the p, and obeying 

M 2 (n) = M 2 p (l + an) , (1.2) 

reproduces reasonably well the parton-model logarithm of the two-point vector corre- 
lation function if a = 2 and the p decay constant is taken at its experimental value. 
This was an exciting result because the Veneziano model ||, which can be considered 
as the model paradigm of Regge theory, predicted exactly this value. We will later on 
see that, if the distribution of eq. ( |1.2|) is assumed, the value a = 2 is actually an exact 
algebraic consequence of the OPE, independent of the value of the p decay constant. 
However, in that case eq. ( fL.2|) also leads to a negative value for the gluon condensate, 
a result which is clearly disfavored phenomenologically U . In our approach we will not 
take the p to be part of the equally spaced tower. 

Also Geshkenbein considered JL0J the problem of an infinite set of resonances in 



the vector channel and how this would match onto the OPE; this time with a gen- 
eral function My(n). However, nobody seems to have considered an infinite set of 
resonances in the vector and axial-vector correlation functions on the same footing, 



together with their properties under chiral symmetry [11|. In addition, the experimen- 
tal knowledge of the vector spectrum has improved considerably recently, so that an 
updated comparison is interesting. 

In this paper, we will study the phenomenological consequences of our ansatz, in 
which we take equally spaced towers in both the vector and axial- vector channels with 
separate resonances for the p and the pion, to leading order in N c , and in the chiral limit. 
Furthermore, although we will refer explicitly to the p and ai mesons for simplicity, it 
should be clear that our results also apply in the U(3) x U(3) chiral limit. 

^n principle one could also play with the n dependence of the resonances' decay constants (see 
below) . 
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The rest of the paper is organized as follows. Section 2 is devoted to the presentation 
of our model. Section 3 contains the phenomenological consequences of this model. In 
section 4 we make a comparison with the Lowest-Meson Dominance approximation of 



ref. [121, and we conclude in section 5. 



2 The model 

Let us start by defining the vector two-point function as 

TL%A=iJ d A xe^{J^ A {x)J% A {$)) = (q,q v - g, u q 2 ) Yl VA (q 2 ) , (2.3) 

with Jy(x) = d(x) r y' J 'u(x) and J%{x) — d(x)^rfu(x). Both functions YIv,a{Q 2 = ~Q 2 ) 
satisfy the dispersion relation (up to one subtraction) 

2 f°° dt 1 



KvAQl = / 7-7^ -Imn^(t) . (2.4) 

JO t + Q z 7T 



Following the discussion in the introduction we will assume that 



oo 



1 lmU v (t) = 2F 2 5(t - M 2 ) + 2 V F 2 5(t - M 2 (n)) , i2.Vi 



n=0 



1 oo 

-ImIT A (t) = 2F 2 S(t) + 2 £ F\5{t - M\{n)) . (2.6) 
* n =o 

In these expressions, F p is the electromagnetic decay constant of the p, and F n = 
93 MeV is the pion decay constant. F v and F A = F ai are the corresponding constants 
for the towers of vector and axial- vector mesons. Note that, along with the ansatz of 
eq. C P- - ID f° r My A (n), we take Fy A both to be independent of n. 

Since QCDoo has not been solved, it is unclear what values to take for all these 
parameters even in the large- iV c limit. We estimate that 1/N C subleading corrections 
are roughly of the order of Y/M ~ 20 — 30% so that decay constants and masses are 
uncertain by this amount. In addition, we will give results for both F n ~ 87 MeV and 
F^ ~ 93 MeV as a rough estimate of how much our results change as a consequence of 
chiral corrections. 

The analysis of this ansatz is just a straightforward application of the Euler- 
Maclaurin summation formula. For instance, in the vector case one obtains 

y *9 = r +l dn 3 + 



Q 2 + nA 2 v + m 2 v Jo Q 2 + nk 2 v + m v 

+ - S + J F v A v + (2 7) 

2 Qi + ml^ 12 (Q 2 + m 2 y ) 2 ' K ' } 

where N v ^> 1, with the condition that A 2 utoff = (Ny+l)ky+my is a momentum cutoff 
which regularizes the theory and is taken equal to the momentum cutoff in the axial 
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channel in order to avoid spurious chiral symmetry breaking in the deep ultraviolet. 
Performing the integral and expanding in inverse powers of Q 2 , one finds 



o rp2 A 2 -i 

TT (D 2 \ — V lr,<r cutoff | 1 



+ 



+ 



1 

Q 4 
l 

0° 



^FM + i^A 2 , 



2F 2 M 4 - -FyAy 




A 6 



3 m v 1 m 2 / 
2 + 2 A 2 " 



(2.f 



+ o(i/q 8 



and in the axial channel 

n A (Q 2 ) = 
+ 



2F\ 
1 

l 



log 



A 2 

^cutoff 



Q 2 



1 

+ g 2 



2F 2 - 2F 2 



F 2 A 
2 



Ai + 6, 



Ai 



A^ 



3 m 4 



2 A 4 + 2 A^ 



(2.9) 



Of course the expansions of Uv,a{Q 2 ) in eqs. 



+ 0(1/Q S ) . 
2~P|) have to match their respective 



OPE expansions [13 . This leads to the constraints 



2 N r 



Ay 



3 

Ai 



3167T 2 

from the coefficient of the parton-model logarithm in front of the unit operator 



(2.10) 



Fy 



Av 



and Ft 



m A 



(2.11) 



from the absence of a dimension-two operator in the OPE, 



12tt 



-2F 2 M 2 + F 2 A 



F 2 A A\ 



m A 

M 




rn 



(2.12) 



from the contribution of the gluon condensate, and 



28 
~9~ 



44 — 



2f;m; 



-I^a 4 



-FyAy 




m A 

Ai 



3 m v 

v , 

2A 4 V + 2A^ 



+ 



1 m; 



2 A 2 / 



3 m 4 



(2.13) 



from the four-quark condensate. 
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There are two more constraints which we can impose using physical observables 
which depend only on Uy^iQ 2 )- These are the coupling L w of the 0(p 4 ) chiral La- 
grangian ]H] and the electromagnetic pion-mass difference. They are expressed in 
terms of liy a{Q 2 ) as 

1 d 
AdQ 2 

and 



Lio = (Q 2 n LR (Q 2 ) o2=n , (2.14) 



where 



m n+ - mn0 = / d Q 2 Q 2 U LR (Q 2 ) , (2.15) 

8irm n F£ Jo 

U LR (Q 2 ) = \ (ny(Q 2 ) - U A (Q 2 )) . (2.16) 

Given F^, the gluon condensate and the combination Tra s (ipif}) 2 , we can solve for 
the eight parameters F p , M p , F v ,a, mv,A and Ay t A from eqs. ( |2.10| - f2~T3| ), and from the 
solution determine L 10 and m n + — m^o. Thus, the system is overconstrained. 



3 Phenomenological results 

In principle, one could solve the system of equations eqs. ( |2.10| - [2.13"D as indicated at 
the end of the previous section. However, the condensates are relatively poorly known, 
and we prefer to follow a different approach. We will use as inputs the values of 
F n = 93 MeV (we will also consider the value F n = 87 MeV, which is the value in 
the chiral limit 0), M p = 770 MeVQ, and F p ranging from about 120 to 150 MeV, 
the latter being the experimental value. (For larger values the combination 7ra s (ipip) 2 
becomes negative, in contradiction with the phenomenologically most favored value. 
See below.) Choosing a range like this is reasonable, because we are in the large- iV c 
limit, and expect variations of order 20 — 30%. It allows us to eliminate the condensates 
from eqs. ( 2.12 , 2~T3| ), and solve the remaining two equations along with eqs. ( 2.1C . 2~TT| ) 



for F v< a, my,A, and Ay ^, and then calculate the two condensates in terms of these. 

The solutions for Fy t A are shown in figure 1, those for mv,A m figure 2, and the 
condensates are shown in figures 3 and 4, as a function of F p . Below F p w 123 MeV the 
equations do not have a real solution. In the axial channel, rriA = M ai = 1230 MeV 



experimentally, while for Fa = F ai a value of 135 MeV is reported in ref. |T1|. Fy 
is not known experimentally. Phenomenological estimates for the condensates are 
a s (G 2 ) = 0.048 ± 0.030 GeV 4 0, and na s ^) 2 = 9 ± 2 x 10~ 4 GeV 6 0. Looking at 



figures 1 to 4, we see that very good agreement is obtained for values of F p of about 
130 — 135 MeV. For F p = 130 MeV we summarize both our solution along with the 
experimental values in table 1, in which we also give our results for F n = 87 instead of 
93 MeV. 



2 Since one should actually use large- N c values, one could try to vary the p mass as well. However, 
since our choice M p ~ 770 MeV already does a good job, we decided to keep the analysis as simple 
as possible and not to fiddle with M p , while we wait for more data to come in (for instance on the 
spectrum) . 
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Figure 1: 

21] ) as a function of F p in eq. (|2~5l) . All in GeV. 



We may now use the values obtained for mv,A arid Ay.A in order to calculate the 
vector and axial- vector meson spectrum, from eq. ( |1 . 1| ) . The results are shown in table 
2. While experimental confirmation for the mass of the a[ (and higher resonances) is 
needed, table 2 shows a good agreement between the model and the real world. 

Using these results, eq. ( |1 . 1|) can now be rewritten as 



n 



Ml (1 + on) 



0.7 



(3.17) 



This is actually very close to the original proposal of ref. || [7| where the p is not a 
separate resonance from the equally spaced tower but just another member. In order 
to make the comparison more clear, rewrite eq. ( |1.2| ) in terms of the p' mass instead of 
the p. Equation ([L2|) (with a = 2) then reads M v (n) = Mj(l + 2n/3) and a = 2/3, to 
be compared with eq. ( |3.17| ). However, although the spectrum is quite similar, other 
quantities are not. As a matter of fact, note that the result a = 2 in eq. ( |1.2|) (i.e. 
m v /A v = 1/2), can also be derived as a consequence of the lack of dimension-two 
operator in the OPE: the solution a = 2 follows from eq. ( 2.11]) (after setting F p = 0, 
since the p meson belongs to the equally spaced tower of refs. [§, [7|). However, eq. 
( |2.12|) then yields a negative gluon condensate (the gluon condensate is very sensitive 
to variations in meson masses and decay constants). 
From eq. ( |2.14| ), we find for the ES ansatz 



'10 



Ml + \K\ 




n 



3 16vr 2 bg F v 



(3.18) 



For F p = 130 MeV the values we obtain are recorded in table 1. The value of L w varies 
by less than 15% over the range of F p we considered. 

For m 7r + — m n o the expression following from eq. ( |2.15| ) is rather lengthy, and we 



relegate it to an appendix. Again, the values for F p = 130 MeV are presented in 
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Table 1: Results obtained for two different values of F n . M p = 770 MeV is input, and 
we choose F p = 130 MeV. 





F n = 87 MeV 


F w =93 MeV 


Experiment 




F p (MeV) 


130 


130 


153 ±4 P] 




M p (GeV) 


0.77 


0.77 


0.7700 ± 0.0008 


|15| 


F v (MeV) 


140 


130 


no available data 


m v (GeV) 


1.44 


1.42 


1.465 ±0.025 | 


g 


A v (GeV) 


1.21 


1.17 


see Table 2 




F ai = F A (MeV) 


160 


170 


135 ±30 (TB 


! 


M ai = m A (GeV) 


1.27 


1.35 


1.230 ±0.040 




Aa (GeV) 


1.42 


1.50 


see Table 2 




a s (G 2 ) (GeV 4 ) 


0.01 


0.02 


0.048 ± 0.030 


1 


vra s (^) 2 (10~ 4 GeV 6 ) 


7 


9 


(9 ±2) 17] 




171^+ — m w a (MeV) 


4.2 


4.7 


4.5936 ± 0.0005 


PI 


L 10 (fi = M p ) (10- 3 ) 


-5.2 


-5.6 


-5.13 ±0.19 [ 
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Table 2: Vector and axial-vector mass spectra in GeV with the corresponding experi- 
mental values, where available. 







pf 


P" 


ft" 


p"" 


a\ 


a[ 


< 




= 87 MeV 


1.4 


1.9 


2.2 


2.5 


1.3 


1.9 


2.4 


F n 


= 93 MeV 


1.4 


1.8 


2.2 


2.5 


1.3 


2.0 


2.5 


Exp. 


P, 18, IS] 


1.47 


1.70 


2.15 


no data 


1.23 


1.64 [IS], li 


\ [1£] no data 
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Figure 2: 

Resonance masses in eq. Ql.lj) as a function of F p in eq. (|2] 



All in GeV. 



table 1. In the calculation of the pion-mass difference we take M v = 135 MeV. At 
F p = 150 MeV, we find that - = 3.0 resp. 3.2 MeV for F w = 87 resp. 93 MeV. 
It is clear that the ES ansatz favors a lower value of F p of around 130 MeV. Since we 
find no solution for F p < 120 MeV, the solution is quite constrained. 

In principle, our model may be used to predict hitherto unknown values of quantities 
of phenomeno logical interest. Since we have restricted our attention to vector and axial- 
vector two-point functions, there are not many such quantities. One such quantity is 
a K — > tt matrix element of the electro-magnetic penguin operator Qj^\ which can be 
expressed solely in terms of Ulr pOR : 



Q 3 /V) 



(tt+IQ 3 / 2 \K + ){^) 



8ttV. 2 



dQ 2 Q A U LR (Q 2 



(3.19) 



Even though this quantity by itself is of limited phenomenological interest, it is interest- 
ing to compare the value we obtain with that obtained from Lowest-Meson Dominance 
fll"2f (LMD; cf. section 4 below). The difference in numerical values gives an indication 
of the robustness of this sum-rule approach with respect to variations in the ansatz 
taken for the spectral functions. In addition, lattice ]2T| and other phenomenological 



estimates |22| [23[ exist p4|| . Note that this quantity is scale dependent, because the 
integral is logarithmically divergent. This is because terms up to l/Q 4 cancel in 
(cf. eq. ( |2.16| )), but not terms of order 1/Q 6 . We will return to this quantity in the 
next section. 



3 The physical K — > (tttt)i = 2 matrix clement is related to this simpler matrix clement by ChPT. 
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Figure 3: 

The gluon condensate a s {G 2 ) (GeV 4 ) as a function of F p (GeV) in eq. ( |2.5|) . 

4 Comparison with the Lowest- Meson Dominance 
approximation 

It is interesting to compare the ES ansatz which we have been considering thus far with 



the Lowest-Meson Dominance (LMD) ansatz of ref. |12 |. In the ES case, the (lead- 
ing) perturbative values of Uv,a{Q 2 ) for large Q 2 come from the assumed Regge-like 
spectrum. One may also directly incorporate pQCD into an ansatz by approximating 

H-vAQ 2 ) with the 

expression one finds in pQCD which, to leading order in a s , reads 

-lmU v (t)\ pQCD = l^6(t-s ) . (4.20) 

Here a parameter so is introduced which determines the onset of the perturbative 
continuum, above which we expect pQCD to give accurate results. For t > s this 
expression replaces the sums in eqs. fl2.5|J2.6| ). All resonances for t < s are kept; for an 
illustration, see figure 5. In the simplest version of this approach, one keeps one reso- 
nance each in the vector and axial- vector channels. Since, by definition, sq in eq. (|4.20j ) 



marks the onset of the perturbative continuum, it is taken to be the same for both vec- 
tor and axial-vector channels, because chiral symmetry is not broken in perturbation 
theory. This ansatz has recently been shown to yield a successful description of Aleph 
data for the moments of the IL^ function [25|| , it successfully predicts pseudoscalar 
decays [^] and, also recently, it has allowed an analytic calculation of the K° — K° 
mixing parameter, [p7| - 

Qualitatively one can understand the transition between ES and LMD in the follow- 
ing way. Low-energy observables are dominated by the lowest individual resonances 
to the extent that the effect of the continuum is suppressed by the fact that s is 
larger than the mass (squared) of these individual resonances. For the LMD case, 
the continuum ( [4.20|) produces a logarithm just like the integral in eq. (|2.7| ) does, 
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with s = 24ir 2 F 2 /N c = m 2 v (see refs. [[T| g§ for details on the LMD ansatz). Tak- 



ing also into account the leading correction term in eq. ( |2.7| ) to order 1/Q 2 , we find 
s = my (l — \A v /m v ^j. Taking A v /m v smaller (keeping F 2 /A v fixed, cf. eq. ( |2 . 1 0| ) ) 
clearly approaches the perturbative approximation of eq. (^4.20|) , with so — > m v , and 
with the correction terms on the right-hand side of eq. (|2.7| ) becoming smaller (for 
any value of Q 2 ). In this sense, the LMD approach can be seen as following from the 
ES approach in the limit Ay /my — > 0, at least in the vector channel. (In the axial 
channel, a similar argument would apply if we had taken one separate resonance apart 
from the tower in the ES approach. In our case, rriA is relatively small, and \A\jm\ 
is substantially bigger than ^Ay/my.) 

Of course, it is easy to put a s corrections into eq. (|4.20|) in the LMD approach, 
while this is less straightforward in the ES approach (it requires a more complicated 
n dependence of M VA and F VA ). One can also improve the extent to which LMD 
approximates ES by keeping more individual resonances explicit in the LMD ansatz. 
This will move up so or equivalently m v , hence reducing the size of A v /m v . 

Using the same information from the OPE, and the relation F 2 = 2F 2 ||29|| , one 
finds that in the LMD approach so is determined to be so = (AnF^ 2 from the vector 
channel, and then, from the axial- vector channel, that M ai is close to, but smaller than 
y/sQ ||28||.p| This implies that within LMD the a>i resonance has to be kept, and cannot 
be folded into the perturbative continuum, eq. fl4.20|) . Higher resonances are well above 
y/so (whether we take their experimental values, or those predicted from ES, cf. table 
2). We will therefore compare the phenomenology of the ES ansatz, summarized in 
tables 1 and 2, with that of the LMD ansatz with one explicit resonance in each channel 
(in addition to the pion). The LMD results are presented in table 3. All quantities 
displayed in table 3 are predicted in terms of F n with the LMD ansatz []12| pgfl , using 



4 In comparing with ref. p8[, note that F„ = f p M pi etc 
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ImfI v , A 



Figure 5: 

Schematic view of the spectrum (set of vertical lines) together with the parton model, 
eq. (|4.20|) (flat horizontal line). 

the relation F? = 2F%, but (obviously) no predictions are made for excited states in 
both channels in this case. In the comparison with the results in table 1, one should 
keep in mind that the relation F 2 = 2F% was not used with the ES ansatz. However, 
for F n = 87, resp. 93 MeV, F p = 130 MeV satisfies this relation at the 10, resp. 2% 
level. 

Comparing the results of ES and LMD, it is clear that both approaches yield consis- 
tent results, especially if one takes into account that all results are in the large- N c and 
chiral limits, so that deviations of less than 30% are not significant. Possible exceptions 
to this are F ai and 7ra s {ipip) 2 , for which the LMD values are smaller, resp. larger, than 
the ES value. Unfortunately, the experimental values, while closer to the ES values, 
are rather poorly known. Possibly with these two exceptions, the differences between 
tables 1 and 3 give a good indication of the spread one can expect when one combines 
QCD techniques such as large N c , OPE and pQCD with an ansatz for the spectrum in 
order to obtain explicit solutions. 

Finally, we compare the ES prediction for Q^ 2 (/i) with the LMD prediction for the 
same number. First, we show the function Q 4 Hlr{Q 2 ) i n figure 6, for both the ES and 
LMD ansatze. We evaluated the ES value of Q^ 2 (jj) at fi = 2 GeV from the area under 
the curve using \i as a sharp cutoff (choosing F p = 130 MeV and M p = 770 MeV), and 
quote the LMD value from ref. |2"0[] , finding, 

f -0.021 (GeV) 4 , ES, F w = 93 MeV, 
Q 3 7 /2 (fi = 2 GeV) = I -0.019 (GeV) 4 , ES, F n = 87 MeV, (4.21) 
[ -0.024 (GeV) 4 , LMD. 
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Table 3: LMD results obtained for two different values of F 7 





F n = 87 MeV 


F n =93 MeV 


F p (MeV) 


122 


131 


M p (GeV) 


0.77 


0.82 


F ai (MeV) 


87 


93 


M ai (GeV) 


1.08 


1.16 


a s (G 2 } (GeV 4 ) 


0.01 


0.02 


7ra s (^) 2 (l(T 4 GeV 6 ) 


13 


19 


m n + — m n o (MeV) 


5.2 


6.0 


L 10 (n = M p ) (10- 3 ) 


-4.8 


-4.8 



We took again F 2 = IF 2 in the LMD case, in which case the LMD value is independent 
of F.,, (for a fixed M p ). We see that, as in the case of most other quantities, the predicted 
value we obtain for Q^ 2 (^ = 2 GeV) is quite robust to the variation of the ansatz for 
the vector and axial-vector spectral functions we have been considering. The values 
we obtain are fully consistent with each other within the context of the large- N c and 
chiral limits. 




-0.00175 



Figure 6: 

The function Q^Hlr^Q 2 ) (GeV 4 ) for our equal-spacing model (ES) (solid line) and its 
lowest meson dominance version (LMD) (dashed line) as a function of the euclidean 
momentum Q 2 (GeV 2 ). 



5 Conclusion 

In this paper we have considered the vector and axial-vector current two-point func- 
tions in QCD. The combination of taking the large-iV c and chiral limits, the use of 
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asymptotic freedom and the OPE constrain the possible form of these two-point func- 
tions. However, this is not enough to obtain concrete predictions: more knowledge, or 
lacking that, an ansatz for the spectral functions in these two channels is still needed. 

A particular ansatz which has been considered in the literature is the LMD ansatz, 
in which the high-energy part of the spectral functions is calculated in pQCD, while 
one (or a few) explicit resonances are kept at low to medium energy. This approach has 
been quite successful in reproducing known phenomenology, and should therefore be 
useful in estimating as yet unknown hadronic quantities, such as weak matrix elements. 
However, although LMD yields a rather intuitive picture it is not a simple matter to 
give it a fully systematic treatment. Such systematic treatment should, in particular, 
produce a theory of the so-called quark-hadron duality violations [30 1 , a problem which 
is clearly beyond the scope of the present work. Consequently, determinations of 
the systematic errors associated with this approach have to rely at some point on 
phenomenological experience. 

One way to address this issue is to vary the nature of the ansatz within the confines 
of our general picture of the strong interactions, and this is what we have done in 
this paper. Instead of LMD, we have based our ansatz on Regge theory, except for the 
lowest states in each channel (the p meson and the pion). For the precise form of the ES 
ansatz, see eqs. (p~5|j2T^ ) . The hadron phenomenology calculated with this assumption 
is again very good, and agrees well with what is found with LMD (cf. tables 1 and 
3). Possible exceptions perhaps are F ai and na s {ipip) 2 ; experimentally these quantities 
are not that well-known. In general, because of the use of large- N c and chiral limits, 
agreement at the 20-30% level, between ES and LMD, or with experiment, should be 
considered satisfactory. 

The gluon and four-quark condensates are rather sensitive to the values of other 
quantities (masses and decay constants). Within our solution (cf. table 1), we find that 
the p comes out to be very close to lying on the daughter trajectory predicted by Regge 
theory in that M%, ~ M 2 + A v and F p rs Fy. However, requiring the p to be precisely 
on this trajectory leads to a negative gluon condensate. Since large N c (supplemented 
with Regge Theory) and pQCD suggest that asymptotically the tower of resonances 
has to be equally spaced, it makes intuitive sense to model the resonance spectrum 
with an equally spaced tower above 1 GeV, keeping those below 1 GeV separate. The 
ES ansatz goes beyond LMD in that it yields a value for the spacing between the 
resonances (the parameters A V a), and leads to a reasonable resonance pattern in the 
vector and axial- vector channels (cf. table 2). 

We have then used the ES ansatz in order to predict a weak matrix element, Q^ (p), 
which is of interest in the context of non-leptonic kaon decays. We find a value rather 
close to that predicted from LMD, thus lending support to the credibility of this pre- 
diction. 

Finally, we discussed the extent to which the two ansatze are equivalent. We qual- 
itatively argued that in the limit in which we take the spacing between the resonances 
in the ES tower to zero one recovers the LMD approach. However, with LMD, we kept 
both the p and the a x as separate resonances (because one finds that the perturbative 
threshold sq is larger than M% ), while we did take the a± to be the first resonance of 
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the ES tower. 

After all these years we still do not have a real understanding of the large-iV c 
expansion in QCD. Clearly the ES model is not meant to be a substitute for large- N c 
QCD at the theoretical level but, given its phenomeno logical success, it may help in 
our understanding of this limit of QCD. 
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Appendix 



The expression for the electromagnetic pion mass difference, following from the ES 
ansatz follows from eq. ( |2.15| ), is 
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This expression results from calculating the integral in eq. (|2.15 ), which is finite, be- 
cause Hlr(Q 2 ) ~ 4:7ia s (ip4>) 2 /Q 6 for large Q 2 . By partial integration 



/•oo 1 /-oo A 

j o dQ 2 Q 2 T1lr(Q 2 ) = -gjf dQ 2 Q^U LR (Q 



(5.24) 



The sums over n inside -^zTIlr{.Q 2 ) now converge. Next, one would like to make 
use of the fact that terms up to order 1/Q 8 in -^Hlr{Q 2 ) vanish, as follows from 
eqs. ( |2.10| - prT2| ). This can be done by writing 
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(5.26) 



iq \p(p + q)...(p+(£- l)q) (p + (N + l)q) . . . (p + (N + 

with N = oo and £ = 1. This process can be repeated by employing this equation for 
higher values of I until the "left-over" sum on the right-hand side of eq. ( [5.25 ) is of 
order 1/Q 8 , so that the integral over Q 2 in eq. ( 5.24] ) can be done term by term. The 
remaining terms combine with similar terms in the axial channel, and the Q 2 integral 
over these combined terms is finite as well, by virtue of eqs. ( |2.10| - p7L2"|) . One ends up 
with the expression given in eqs. (|5.22|,|5.23D . 
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